Spin effects in transport through non-Fermi liquid quantum dots 
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The current-voltage characteristic of a one dimensional quantum dot connected via tunnel barriers 
to interacting leads is calculated in the region of sequential tunneling. The spin of the electrons 
is taken into account. Non-Fermi liquid correlations implying spin-charge separation are assumed 
to be present in the dot and in the leads. It is found that the energetic distance of the peaks in 
the linear conductance shows a spin-induced parity effect at zero temperature T. The temperature 
dependence of the positions of the peaks depends on the non- Fermi liquid nature of the system. For 
non-symmetric tunnel barriers negative differential conductances are predicted, which are related to 
the participation in the transport of collective states in the quantum dot with larger spins. Without 
spin-charge separation the negative differential conductances do not occur. Taking into account 
spin relaxation destroys the spin-induced conductance features. The possibility of observing in 
experiment the predicted effects are briefly discussed. 

PACS numbers: 73.63.Kv, 71.10.Pm, 72.25.-b 



I. INTRODUCTION 

Spin phenomena in the transport properties of low- 
dimensional quantum systems have become a subject 
of increasing interestiS Several important effects have 
been found when controlling transport of electrons one 
by one in quantum dots. For instance, a parity effect 
in the Coulomb blockade^ of quantum dots with very 
small numbers of electrons, and in Carbon nanotubes^ 
have been detected. The spin blockade effect, especially 
in the non-linear current-voltage characteristic of one- 
dimensional (ID) quantum dots has been predicted££ 
Combining the spin blockade with spin-polarized detec- 
tion, the electron spin in a lateral quantum dot has been 
probed, and spin-related phases were found that have 
been associated with correlations between the electrons^ 

The former parity effect is related to the Pauli principle 
and is quantitatively affected by the contribution of the 
exchange interaction towards the energy of the ground 
state of the electrons occupying the quantum dot. The 
spin blockade effect, which leads to a negative differen- 
tial conductance (NDC) in the current as a function of 
the bias voltage, is due to a combined influence of the 
exchange contribution to the energies of the correlated 
electronic eigenstates and the spin selection rules for the 
transport processes. For example, an excited n-electron 
state of a quantum dot generally can be depopulated via 
two spin channels, namely either by increasing or by de- 
creasing the z-component of the total spin, s, by 1/2. 
However, when the state with the highest total spin be- 
comes populated at a certain bias voltage V, it can be 
depopulated only via processes that decrease s. This can 
lead to a reduction of the total current I(V), when in- 



creasing the bias voltage, thus yielding 81/ dV < 0. 

This phenomenon has been predicted by using a simple 
model for the transport mediated by sequential electron 
tunneling processes through a ID quantum dot contain- 
ing few electrons (n < 4). The electronic eigenstates 
of the latter have been determined numerically exactly 
in the presence of interactions and spini&i The tunnel- 
ing matrix elements of the barriers connecting the states 
of the quantum dot to those in the Fermi-liquid leads 
were assumed to include spin selection rules via Clebsch- 
Gordan coefficients and assuming ad hoc some tunneling 
coupling matrix elements between the correlated quan- 
tum dot states and the states in the leads. Only later, 
the matrix elements have been calculated microscopically 
in order to test the former assumptions £ There are sev- 
eral, in the details somewhat different and intricate mech- 
anisms that produce such a spin-induced reduction of 
the current. Characteristic dependences on the magnetic 
field may be used to distinguish between them. While 
certain signatures of these features have been found in 
a recent experiment done on 2D quantum dotj^S exper- 
imentally well-controlled evidence in quasi-lD quantum 
dots is missing. 

The recent experimental realization of semiconductor- 
based ID quantum wiresiiii2ii£ has opened new perspec- 
tives to systematically investigating the influence of in- 
teractions, spins and impurities on electron transport 
properties. Signatures of spin-charge separation have 
been observed^* and analyzed^ in the tunneling between 
parallel one-dimensional wires. Also carbon nanotubes 
can now be controlled to such a high degree that in- 
vestigations of electronic transport features have become 
possible. 4 i 5 i 16 i 17 The effects of non-Fermi liquid correla- 



2 



tions in ID quantum dots have been analyzed non per- 
turbatively in the coherent tunneling regime^iS 

In the non-linear transport spectra of these devices, 
obtained as the derivative of the current-voltage char- 
acteristics at different gate voltages, a large number of 
low-energy excited states have been founds These can- 
not be understood only in terms of charge excitations 21 
but are theoretically predicted to be related also to the 
spimSS Thus, one can expect that ID spin blockade effects 
could be seen in an experiment. Among the experimen- 
tal realizations of ID quantum dots are electron islands 
between two successive impurities in a ID quantum wire 
(containing interacting electrons) and carbon nanotubes. 
Thus an extension to the theory of the blockade effects 
should include not only the generalization to higher elec- 
tron densities — the previous calculations have been done 
in the limit of low electron density — but treat the inter- 
actions and the spins within the quantum dot and within 
the leads on an equal footing. 

There are perspectives for application of quantum 
structures in spin-electronics, quantum computing and 
communication,^ Previous works focused on spin trans- 
port in (2D) quantum dots connected to non-interacting 
leads in the presence of a magnetic field, 24 includ- 
ing an oscillating magnetic electron spin resonance 
component^ Spin transport in circuits with ferromag- 
netic elements and Luttinger- liquid interactionSSiSL 2 ^ has 
been considered. In view of the applications, the theory 
of the spin control of electron transport in the presence 
of correlations is very important since in nanoscale de- 
vices the latter can be very important. In our previous 
work^ we have found that at T — 0, spin polarization 
effects are robust against the correlations. They can even 
be enhanced, if the polarization is not complete. 

In the present paper, we extend our previous calcula- 
tions to include temperature effects, asymmetry in the 
tunnel barriers and the effect of spin relaxation. The 
former are treated microscopically within the sequential 
tunneling approximation using the microscopic model of 
a ID quantum dot described by a Luttinger liquid of fi- 
nite length connected via tunnel barriers to semi-infinite 
interacting leads. Spin relaxation is treated phenomeno- 
logically by including tunnel rates corresponding to spin- 
flip transitions into the master equation for the dynamics 
of the probability distribution. Our main results are 

1. We deduce that the temperature dependence of the 
linear conductance peaks reflects the non-Fermi liq- 
uid correlation in the leads. 

2. In non-linear transport, we find that many-particle 
states with higher spins in the quantum dot act as 
traps such that the electric current decreases with 
increasing bias voltage. The physical origin of this 
is spin-charge separation. 

3. Spin relaxation leads to a destruction of the neg- 
ative differential conductance peaks. This con- 
firms that they are related to the trap properties 



of higher-spin states. A typical "phase diagram" of 
the crossovers between negative and positive differ- 
ential conductance peaks is calculated. 

We discuss the quantitative conditions for observing the 
predicted effects in an experiment. 

In the next section, we introduce the model for the 
Luttinger liquid quantum dot. Section 3 is used to dis- 
cuss the energy scales of the system. In the section 4 we 
describe the approach for calculating the transport prop- 
erties. Section 5 is devoted to the linear and section 6 
to the non-linear regime. In section 7 we compare with 
other recent approaches, discuss possible experimental 
realizations and draw some conclusions. 



II. TOMONAGA-LUTTINGER MODEL FOR A 
ID QUANTUM DOT WITH SPIN 

We consider the ID system shown schematically in 
Fig-ffl The region |x| < a/2 represents a finite quantum 
wire that plays the role of the ID quantum dot. This is 
assumed to be connected to reservoirs with electrochem- 
ical potentials /i>, = ±eV/2 (A = R, L) via semi-infinite 
ID systems between — L < x < —a/2 and a/2 < x < L 
(L — > oo ) that represent left (L) and right (R) leads, re- 
spectively. The electrochemical potentials are controlled 
by the bias voltage V. Tunneling barriers at x = ±a/2 
connect the quantum dot with the leads. The three 
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FIG. 1: Schematic representation of the Luttinger liquid 
quantum dot (|a;| < a/2) connected via tunneling barriers 
to semi-infinite wires. The charge (g p ) and spin (g a ) inter- 
action parameters in the leads and the dot are assumed to 
be different. The bias voltage V controls the chemical poten- 
tial difference of the leads. The gate voltage V g controls the 
chemical potential quantum dot (cf. Fig. |2J. 

regions of the quantum wire are assumed as interact- 
ing Luttinger liquids (LL) 29 i 30 ' 31 i 32 with possibly differ- 
ent interaction constants for charges (<? p ) and spins (g a ) 
(Fig.nj. 33 34 In addition to the bias voltage V we assume 
that the chemical potential in the dot region is tuned via 
a gate voltage V s . The total Hamiltonian is 

H = H + H t +H c . (1) 

Here, H a = + + describes the three un- 
coupled LL, H t the tunneling term, and H c the coupling 
with the external electric circuit. 
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The ID quantum dot (|x| < a/2) is treated with 
open boundary conditionSfi&Si i/ji d \±a/2) = 0, with 

the fermion operators ipi (x) corresponding to the z- 
component of the spin s = ±1 (unit h/2). 

The Luttinger liquid theory describes the energetically 
low-lying excitations of ID interacting electrons around 
the Fermi energy Ep. It relies on the linearization of the 
dispersion relation around the Fermi level giving rise to 
a bare, constant level spacing eo = tmvp/a. As a nec- 
essary condition for this assumption, one needs to fulfill 
the condition Ep ^> e . For a system of finite length, this 
implies a constraint between the Fermi wave number and 
the momentum discretization, kp 3> ir/a. There is also 
a constraint with respect to temperature induced by the 
finite size of the dot. In order to resolve the energy levels 
of the correlated states the temperature should be lower 
than the level spacing, k^T <C Eq. 

Near the Fermi point kp the fermion operators are de- 
composed into fields that propagate to the right (r) and 
to the left (1) 



ik F x / (d). 



(2) 



Because of the boundary conditions these fields are not 
independent. The bosonization description is done using 
only tjjs^T (x) . The system is diagonalized in the presence 
of interactions. The corresponding Hamiltonian is2^& 
(units such that H = 1) 



H o d) = E E^(<^ t (<zM<?) + £ 



v=p,a ij>0 



U P -2 

—n 

9p 



— s 2 

9a . 



(3) 

Here, v'(q), v{q) are the boson operators of the collective 
charge (y ~ p) and spin [y = a) density waves (CDW and 
SDW). Due to the boundary condition the wave number 
is quantized, q = irm/a (m integer > 1). The energy 
spectra an 



The zero mode operators h and s represent the excess 
number of charges, and of the z-component of the total 
spin with respect to their average values on the ground 
state. The latter correspond to hq — 2kpa/n — 1 and 
so = in the absence of a magnetic field. The eigenval- 
ues of the zero-mode operators are integers n, s with the 
constraint, due to the boundary conditions, n + s = even. 
The zero-mode energy contributions in (0 represent the 
energy needed for changing the total charge and spin with 
respect to the ground state. 

For the leads we assume a LL with open boundary 
conditions at the tunnel barriers. At L = ±oo they 
are assumed to be connected to reservoirs with different 
electrochemical potentials. The difference of the latter is 
proportional to the bias voltage. In state-of-the-art semi- 
conductor quantum wires (and also in carbon nanotube 
systems) the exchange is only a very small correction as 
compared to the Coulomb part of the interaction^? which 
in any case is expected to be considerably weaker that in 
the region of the quantum quantum dot. The exchange 
interaction in the leads is not expected to influence the re- 
sults described below on the qualitative level. Especially, 
it is not expected to influence the positions of the con- 
ductance peaks that are determined by the interactions 
in the region of the quantum dot. At most, it can give a 
small quantitative correction of the dependence on tem- 
perature of the conductance peaks which is determined 
by the global properties of the system as has been dis- 
cussed earlier^* Therefore, we neglect in this paper the 
exchange interaction in the leads completely for the sake 
of simplicity. The Hamiltonian describing the excitations 
in the leads is thenSi 



H, 



{ X) = E 

v=p,a k>0 
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with the interaction parameters 
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where 



Vn = 



V(0) 
2nvp 



V = 



V{2kp) 
2izvp 



(4) 



(5) 



(6) 



are proportional to the forward (q — > 0), and part of 
the backward (q — > 2kp) contributions of the electron in- 
teraction, respectively. The quantity V(q) denotes the 
Fourier transform of the interaction potential inside the 
dot. The backward term corresponds to an exchange in- 
teraction. The parameters fulfill V cx < 1, necessary in 
order to have a bounded Hamiltonian, and Vq > V CK . 
This implies < g p < 1 (repulsive charge-charge interac- 
tion) and g a > 1. 



with vp the Fermi velocity in the leads, v x {k), v x (k), the 
boson operators of collective CDW, SDW and h\, s\ the 
zero mode operators for the excess of charge and spin 
with respect to their average values in the leads. The 
energy spectra are (fc = Trm/L, m integer > 1) 

njk) = —k, n„(k) = v F k 

9a 

with the charge interaction parameter go 
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2U° 

TTVp 
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< 1 



determined by the average interaction Uo in the leads. 

The coupling between the leads and the dot is de- 
scribed by tunnel barriers at x-^ — —a/2 and = a/2 



E 

s=±l A=L,R 



E M^a)^^) 



h.c. 



(7) 
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with the right moving fermion operators normalized to 
the shortest wavelength, and i L ,R the transmission am- 
plitudes of the barriers^ 

The operator of the external bias and gate voltages 
that allows to electrically control current transport and 
charge density in the dot are written in terms of the op- 
erators h of the excess charges (Fig.|2J 38 



eV 



SC v + ^v. 



h. (8) 



Here, V and V g are the bias and the gate voltage, Cs = 
Cl + Cr . + Cg is the total capacitance with C l ,Cr, and C g 
the capacitances of the leads and the gate, where 8C = 
Cl — Cr. 
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FIG. 2: Equivalent circuit for the quantum dot with asym- 
metric tunnel barriers. Left and right tunnel junctions are 
parameterized by the capacitances Cl,r and the resistances 
Rl,r = LUc/ite 2 ^^, (cf. JlH ), w c cutoff energy of the leads, 
th,R transmission amplitudes of the barriers, C g gate capaci- 
tance, Vg gate voltage, V bias voltage. 



III. ENERGY SCALES 

The states of the isolated quantum dot are 
|n, s, , {'g }) where n and s are the excess numbers 
of charges and spins with respect to the ground state with 
charge uq and spin sq, and {i^} the occupation numbers 
of the CDW = p) and SDW (y = a) at different q. 
Some examples of these states are schematically shown in 
Fig. El Figure^ represents the state |1, 1, {0}, {0}), cor- 
responding to one excess charge n = 1 with spin s = 1. 
Figure^ shows the state |2, 0, {0}, {0}) with two extra 
electrons and s = 0. A higher spin state (s = 2) is shown 
in Fig. Excited states with unchanged n and s con- 
tain CDW and/or SDW. Figure OJi shows an example 
with the smallest possible wave- number qo = it /a and 
smallest excitation energy for n = 2 and s = 0. The 
operator 

^ ( qo ) |2, 0, {0}, {0}> = |2, 0, {0 . . .} , {1, 0, . . . , 0, . . .}) . 

creates a linear superposition of the state (cf. Fig. JSJi)) 
and the one with inverted spins. In the presence of inter- 
action a charge excitation with the smallest wave number 
would be associated with a much larger excitation energy. 
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c) 



FIG. 3: Schematic examples of excited states with respect to 
a background charge no = (2ak F /ir) — 1 and spin so = 0. (a) 
n = 1 and s = 1, Ui = W(l, 1,0,0) (see Eq. ©, (b) n = 2 
and s = 0, U2 = W(2, 0, 0, 0), (c) excited state with s — 2 for 
n = 2, (d) spin excitation with respect to the state (b). 



Excitations with higher wave numbers induced by oper- 
ators v'(q > qo), as well as multiple excitations [i^(qr)] r 
(r > 1) are also possible. 

Using the Hamiltonian © we write the total energy of 
the particular configuration |n,s,{^},{^}) as 



E E 

U (n, s, l p , l a ) = -fn 2 + -fs 2 + l p e p + l a e 



(9) 



The first two terms represent the contributions of charge 
and spin additions. The second two terms correspond 
to the CDW and SDW. Since the excitation spectra are 
linear in the wave number, the energies of the collective 
modes depend on the total numbers of excitation quanta 
l v = J2 q Qlq onr y y i a the discrete excitation energies e v — 
nv v /a. From the microscopic theory one gets (y = p, a) 



EL 



2ag u 



£0 1 + v a 
2 9l 



(10) 



where Eq = irv-p/a the constant level spacing in the non- 
interacting case. These energies are different from zero 
even without interaction due to the discrete nature of the 
energy levels inside the dot and the Pauli principle. 

Despite the microscopic model provides quantitative 
estimates for the charge and spin addition energies, sev- 
eral influences that occur in experimental setups are here 
neglected. For instance, the coupling with the gates, 
long range interaction effects, and the screening due to 
the nearby 2DEG present in cleaved edge overgrowth 
systemsi&Sfi affect E p and can cause strong deviations 
from the simple expression HlDf) , 2 - More robust is the 
spin addition energy which is influenced at most by the 
exchange part of the interaction. Therefore, we treat in 
the following E as a free parameter with E p ^> E a . 

Phenomcnologically, one identifies E p with the total 
electrostatic energy (Fig. 0) E p = e 2 /C s of the circuit 
model. With this one includes the term enV g C g /Cs com- 
ing from © into the zero mode of the charge sector in 

$b,E p (n- n g ) 2 / 2 with n g = VgCg/e- 

For the spin addition energy we continue using the mi- 
croscopic expression (|1(JH 



E„ 



(1 - V cx ) 
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Due to the exchange interaction this energy is decreased 
as compared to the non-interacting case. The energy 
quanta of the collective excitations can be expressed mi- 
croscopically as 



£o- 



1 



V 

r cx 



9u 



(11) 



The presence of the exchange interaction renormalizes 
simultaneously the energy E a of the total spin configu- 
ration and the energy e a of the spin waves 



e* = e a ^l-V 2 x = 2E 



1 + V C , 
l-V„ 



The plasmon energy e p is, on the other hand, affected by 
the Coulomb interaction with in general e p > e a . The en- 
ergetic difference between the CDW and SDW indicates 
the presence of the spin-charge separation in a LL. 

The above discussion implies a hierarchy of energy 
scales characteristic for the model considered 



2E <J < e a < e < Sp < E. 
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This suggests to use E a as the natural energy scale. 



Without interaction 2E„ 



e . 



IV. TUNNELING 



A. Transition rates 



For high tunnel barriers and not too low temperatures, 
the dominant processes that contribute to the electron 
transport are sequential transfers of single electrons with 
spin up or down through the two barriers. In this case, H t 
is treated as a perturbation. Higher order coherent pro- 
cesses can be safely neglected as long as k^T ^> SE, with 
5E the level broadening of virtual states proportional to 
the tunneling ratesi^ We consider processes that change 
the state of the dot from an initial |i) to a final state |f), 

|i> = K s U {l p qi }, {Z£}>, |f) = \n u at, {l p q{ }, {!£}) . 

Because of the sequential nature of the tunnel processes 
these states have to be consistent with the charge and 
spin selection rules 



An = rif — rij = ±1 , As = Sf — Sj = ±1 . 



,(A) 



(12) 



For obtaining explicitly the tunneling rates r[., WOO 
the left (A = L) and right barriers (A = R) we sum over 
all possible final lead states and perform a thermal av- 
erage over the initial states with the chemical potentials 
/iA- From (JJJ one finds 



r (A) 
1 |i>Hf) 



-f2,„A 



dr e lAUT e~ 



(13) 



The matrix element (p^ is 



^d= (i|VS(*A.0)|f> 



The energy difference 

AU=U{\T))-U{\i))- An 



SC ' 



(14) 



is obtained from U defined in @. The signs ± refer to 
the left (+) or right (— ) barriers, the factor exp [— W^t)] 
results from the trace over the lead excitations^ 



leads 



(16) 



The thermal average is performed with respect to the 
decoupled Hamiltonians Hq X K This factor turns out to 
be independent of the spin of the tunneling electron and 
the barrier variables. The contribution of the quantum 
dot is contained in AU and in the matrix elements tp^. 

The tunneling rates rj^^ifi were recently used by Kim 
et alJSL in order to study the out-of-equilibrium dynamics 
in the presence of an external voltage. In this work it has 
been assumed that the collective modes are stable and do 
not relax during the tunneling processes. The probability 
distribution was obtained self-consistcntly using a gener- 
alized master equation. 

The opposite situation has been addressed by Braggio 
et al. 22 who considered the possibility of fast relaxation 
of the collective modes induced by extra processes^ that 
are not included in the diagonal Hamiltonians , such 
as spin-orbit interaction and phonon coupling. In this 
case, the dynamical variables consist of the total charge 
and spin only, while the collective modes are traced out 
with a thermal average over the initial states and a sum- 
mation over the final states. In the following, we use this 
second model. 

As already pointed ou"b^&, the two approaches do not 
give qualitatively different results regarding the current- 
voltage characteristic of the system. However, in a very 
recent work^ the presence of stable plasmons in a spin- 
less ID quantum dot has been shown to affect dramat- 
ically its shot noise properties. On the contrary it was 
also predicted a crossover towards the results of a model 
with fully relaxed bosonic modest, in the presence of 
a phenomenological relaxation rate for the plasmons. 
In particular the two models give the same results if 
7p/ro <^ 1; where 7 P is a phenomenological relaxation 



(L)" 



rate of the plasmonic modes and f q 1 = T 
with g p = g a = 1 - see Eq. 1261 - is a characteristic tunnel- 
ing rate. In a system with spin, we can expect a similar 
result if the phenomenological relaxation rate of the spin 
density waves j s , satisfies 7 s /ro > 1. For semiconductor- 
based ID quantum dotsSi one can estimate Tq « 10 11 
s _1 by using a tunnel resistance i?R, w 100 h/e 2 and 
a w 0.2 /im - see section IY III . It is difficult to evaluate 
microscopically the relaxation rates 7 P , S - However, one 



6 



can estimate^ using a level spacing of about 1 meV, 
7 S » 8 • 10 12 s _1 . It can be expected that j p is even 
larger. In the following we assume 7 P . s /ro > 1 an d thus 
complete relaxation of the bosonic modes. 
The reduced rates are then given by 



-i(A) 

\m,si)— >\m,sf) 



(A) 



10 



(17) 



with P({^;}) the thermal probability distribution with 
respect to Hq' . Performing the sums one finds 



r 



(A) 

| Hi, sj)— >\nf ,Sf) 



dr e ^UT e - Wl (T) e -w d (T) 



(18) 



where 
AU 



E„ 



[1 + 2 (nj 



n g ) An\ 



2 



[1 + 2s;As] 



where the rightmost process (=$ arrow) is associated with 
a fast time scale and the intermediate state contains col- 
lective excitations. In this approximation, only the en- 
ergy of the collective excitation is detectable. In the fol- 
lowing we will use the notation 



^) lp: let) 



(24) 



to label an excited state when it is involved as a final 
state in a tunneling process. 

Let us now investigate in more detail the energy de- 
pendence of the rates. For simplicity, we do not specify 
the initial and final spin and charge states. Exploiting 
the discrete nature of the dot spectral density we rewrite 
(|18fl as a function of the energy difference E = AU 

(E) = lf> % *h, 7 (E - l p e p - l a e a ) , (25) 



where 



-At 



SC , 



is the energy difference associated to the particular pro- 
cess, the signs ± refers to the left (+) or right (— ) bar- 
riers. The kernel e~ Wd ^ represents the thermal aver- 
age over the initial spin and charge collective modes in 
the dot. Using the bosonization method one obtains^ 
(/T 1 = k B T) 



— — — < coth I J [1 — cos {lot) 



W 14 (t) = / d^ 



-i sin (wt) 



(20) 



,(A) 



with 



e 2 r(l + a) 
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(26) 



(27) 



and G\ = R^ 1 = ne 2 t\/uj 2 the intrinsic conductances 
of the barriers. The function 7(2:) is determined by the 
leads^S 



7(s) = l. e ^ 



At T = we have 



2 ff 2tt; 



2tt 

#4 



^ ■ (28) 



Here, the spectral densities of the leads and the quantum 
dot are, 



J\(lu) = —e 
9 



-Lj/u 



and 



j d ( W ) = w e ^E^ 

^e{p,cr} y " m=l 



(21) 



(22) 



The cutoff lo c defines the highest energy in the model and 
.r 1 = (l+% 1 )/2. 

The structure of the spectral density of the dot, J d (ui), 
indicates that, though they are infinitely quickly relaxing, 
the CDW and the SDW still contribute to the tunneling 
dynamics. Even if the relaxation prevents the collective 
excitations to be initial states for the tunneling, it is still 
possible to reach an excited state as a "final" state with 
a given energy. A typical process is 

\n h Si ) -> |n f , at, {l p qi }, {%}) =J |n f , s f ) (23) 



6{x). 



(29) 



The weights ai v are due to the dot contribution. At T 
0, they are 



r(i/2 ff „ + u 
r(i/2^)/j 



0(O 



(30) 



For finite temperature, the weights have to be numeri- 
cally determined. In (|28H - 130|I . T(z) is the Euler-gamma 
function. Figure 01 shows the effect of the charge-spin 
separation in the structure of (E) /T ( Q X) at T = 0. In 
the absence of spin-charge separation (g p = g a = 1) all of 
the energies are degenerate, e p — e a = Sq. It is then easy 
to write down explicitly the contributions to the sum i12">l) 

al9(x) + 2a a 1 9(x - e ) + {a\ + 2a a 2 ) 8(x - 2e ) 
+ 2 (a a 3 + a^) 0(x — 3e ) + • • • 

with do = 1, ai = 1/2, a 2 — 3/8, 03 = 5/16, . . .. Each 
of the coefficients of the theta-functions in the latter ex- 
pression sums up to one (Fig.^J, so one always has unit 
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E 



FIG. 4: Schematic plot of the transition rate T^ X \E), in units 
of Tq A ' , as a function of the tunneling energy E, for go — 1 
and T = 0. Full line: no interactions in the dot, g p = g a = 1; 
dashed line: with interactions, g a > 1 and p p < 1. Excitation 
energies are indicated at the energy-axis. 



steps centered at integer values of £o- If 9a > 1 and 
g p < 1, charge and spin modes are energetically split, 
e a < Eq < e p . Then, the first three terms of Q25|l are 



9{x) 



1 

2^ 



-8{x - Ea-) + ^—0{x - E p ) 



This implies non-integer quantized steps at different en- 
ergies with heights depending on the spin and charge in- 
teractions (Fig. 21 dashed line). For simplicity, in Fig. 0] 
T = and go = 1 have been assumed. Finite tempera- 
tures and go < 1 do not drastically change these results 
but smoothen the jumps at the positions given by the 
dot parameters. 



B. The master equation 

Having assumed fast relaxation of the collective exci- 
tations in the quantum dot necessarily implies that the 
stationary states of the dot are fully characterized by 
the variables n and s. One can then define an occupa- 
tion probability P{n, s) for a state \n, s) that satisfies the 
master equation 

dtP n ,s(t) = J2n'=n±l [Pn',s'('t)^\n',s'}^\ri,s) 
s'=s±l 

} , (31) 

where V = Ylx^ ■ I n the stationary limit, the l.h.s. of 
(|31|l is zero such that one has to solve a homogeneous 
system of linear equations. The solution must be nor- 
malized, P ntS = 1. The stationary current is 



s)^\n+l,s+q) 



-r 



(R) 

\n,s) — >\n — l,s-\-q) 



(32) 



If both eV, ksT < E p at most two charge states en- 
ter the dynamics, and the current is given by tunneling 



events corresponding to transitions |n, s) — > \n + 1, s') — > 
\n,s"). Below, we will consider this regime. 

In the above discussion, it is assumed that states with 
spin \s\ > 1 are stable. Spin-flip processes, however, 
can induce a relaxation of the total spin. These can be 
due to different mechanisms such as magnetic scattering 
and spin-orbit interaction. Another possible source of 
spin-flip processes can be cotunncling. The quantitative 
microscopical evaluation of the corresponding relaxation 
rates is not easy and depends on the dimensionality of 
the system and on the electronic correlations. In GaAs- 
based quantum dots, it seems that spin-flip processes 
have much smaller rates^i, four orders of magnitude or 
even smaller than the non-spin-flip processes involved in 
the relaxation of the bosonic modes. In ID systems, the 
non- Fermi liquid nature of the interaction could lead to a 
non-trivial energy dependence of these relaxation rates. 

It is primarily not the aim of the present paper to study 
these processes. In order to get first insight into the sta- 
bility of the negative differential conductance features it 
is sufficient to introduce a phenomenological spin-flip re- 
laxation rate T^^, , with s' = s ± 2, and modify accord- 
ingly Eq. (13 1 11 - The results of this modification will be 
discussed in section lYl Dl 



V. LINEAR TRANSPORT 

In this section we specialize (|32H to the linear regime 
(V — > 0). We assume symmetric barriers, = £r, and 
Cl = Cr, hence we drop for convenience the barrier in- 
dex in the expression for the rates. The electrochemical 
potential of the quantum dot is defined as 

Md(«. s n) =U{n + l,s„ +1 ,0,0) -U(n,s n ,0,Q) (33) 

where s n , s„+i are the spins of the ground state for n, 
n + 1 electrons, respectively. From @ one finds 

Md (n, S J =E p (n+±-n g ^+ (-1)"^ . 

In the following, we evaluate the linear differential con- 
ductance G = dI/dV\ v=Q . 



A. Low temperature 

At low temperatures, k B T <§C E al the master equation 
is solved in the subspace defined by the ground states for 
ii and n+I electrons. The conductance is 



G(0 



(3e 2 



me 



V8cosh{[/3£ +(-!)" log 2]/2} 



(34) 



where T(£) is the rate in and £ = E p [n g — n g es (n, 0)] 
corresponds to the deviation from the resonance peak 
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position n^ es (n,0) at T — 0. The latter is determined by 
the condition /i^ = 

< s (n,0) = n+i + (-l)»^-. (35) 

The position of the linear conductance peaks is affected 
by the spin that leads to an even-odd effect in the dis- 
tances 6(n <-> n + 1) between conductance peaks 

5(n^n + l) = l + (-l)" +1 — . 

E p 

This even-odd effect was recently observed in carbon nan- 
otube experiments^ 

For T / the peaks in the conductance no longer 
occur at the zero temperature positions. Instead, they 
shift linearly with temperature with slopes that depend 
on the interaction in the leads (Fig. |5j 



>,T)=n~(M) + (-!)"+ V(flo) 



E„ 



Here, <p is obtained from the implicit relation 

"log 2 + </>(#()) 



Im 



, 1 , Xffo) 
1 2g 2tt 



7T 

+ — tanh 



(36) 



(37) 



with ip(z) the digamma-function. Without interaction 
in the leads (g = 1) one recovers the well know result 
0(1) = log 2/2M 




FIG. 5: The prefactor <j> = 2ij> (g ) / log 2 in i'Mii as a function 
of the interaction in the leads. Inset: low-temperature shift 
of a conductance peak for g = 1.0 (solid), g = 0.8 (dashed), 
g — 0.6 (dotted), g = 0.4 (dash-dotted) (temperature units 



For k B T -C e v , since one cannot excite charge and spin 
density waves in the dot, the temperature dependence 
of the conductance maximum is determined only by the 
interactions in the leads 



G max (T) oc T c 



(38) 



B. High temperature 

For E a , e a < k^T the peak position deviates from the 
linear behavior and approaches, for k^T k, E a , the spin- 
less value n + 1/2 (Fig. For k^T ^> E a , one obtains 



n~(0) 




n + 



FIG. 6: Position of the conductance peak for even n, as a 
function of the temperature T (units E a /kB) for g = 0.8, 
g p = 0.3 and g a = 1.0. At high temperatures the peak po- 
sition approaches the spinless value n + 1/2. Inset: double 
logarithmic plot of the conductance maximum G max (T) in 
arbitrary units. Crosses: analytic result from 14011 . 



an analytic expression for the conductance. We factorize 
the occupation probability, P n>s = p(n)^>(s) with \f(s) 
the thermal occupation probability 



*(«) = 



z~ 



-0E„s 2 /2 
-l3E a (s 2 -l)/2 



s even 
s odd 



(39) 



The prefactors Z c = Z Q = ^/ir/2f]E r7 are determined by 
the normalization conditions J2 S ^(2s) = \^(2s+l) = 
1. Using the factorized form for P„ )S we sum over the spin 
states in the master equation obtaining p(n)T n -^n+i = 
p(n + l)r n _|_i_>„ (n even) with 



f n— m+i — / d 

LI «/ 



- s * < 1 |n,2s)— >|n+l,2s+p) 



p=±l ' 



and 



|n+l,2s+p>— >|n,2s> 



The hnal result for the differential conductance is 
/3e 2 r(( + E a /2)T(-<;-Ej2) 



G(0 



2 T(C + E a /2)+T(-C~E a /2) 



(40) 



with C = £ p (n - n g + 1/2). When /c B T > £ CT the con- 
ductance peak as extracted from l|40l) is centered around 
C = which implies n T g S {n,T) = n + 1/2, the spinless 
position. The conductance peak maximum however fol- 
lows different power laws depending on the temperature 
range considered 



with a given in (|27|l . 



G max (T) oc 



j>a-l+(2g a ) 1 
T a-l + (2g a )- 1 + (2g P Y 



e a < k B T < e p 
e„ « k B T < E p 
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Thus, for a well developed energetic separation between 
charge and spin it should be in principle possible to detect 
three different power laws in the linear conductance peak 
maximum value, over the whole temperature range. 

VI. NONLINEAR TRANSPORT 

In the following, we discuss the transport in the non- 
linear regime. We assume n to be even, then we la- 



-1 G 2 




FIG. 7: The differential conductance G (arbitrary units) as a 
function of bias voltage V (units E a /e) and number of gate 
charges n s . System parameters are A = G L /G R = 50, feT = 
10~ 2 E a , C L = 5C* R , Cg/CR = 0.01 g p = 0.63, g a = 1.15, 
<7o = 1) E p = 5E a . Top: gray scale. 

bel the states \n,s) using the spin variable s alone such 
Iu . , , — * r.^, .,. The master equation for the sta- 
tionary probabilities P n . s is solved numerically. From this 
the current-voltage characteristic I(V,n g ) and the dif- 
ferential conductance G(V, n g ) = dl/dV are calculated. 
Figure [7| gives an overview of the behavior of the differen- 
tial conductance as a function of bias voltage V and gate 
voltage induced number of charges n g . Most pronounced 
are the parity effect for V — > and the rich structure due 
to the excited states of the quantum dot. This can be 
understood in detail by considering all the possible tran- 
sitions and taking into account the correlation induced 
features of the tunneling rates. Most important are the 
negative differential conductance peaks indicated by the 
white lines. They are closely related to the non- Fermi 
liquid properties of the model. This is discussed below 
in more detail. 



A. Symmetric barriers 

Figure |H1 shows the differential conductance G(V,n g ) 
for symmetric barriers as a function of bias and gate 
voltages. The conductance is always positive. Each of 
the lines corresponds to a peak in the differential con- 
ductance and is related to a transition between states in 
the dot for n and n + 1. 
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FIG. 8: The differential conductance G (arbitrary units) as a 
function of the bias voltage (units E a /e) and n g , for symmet- 
ric barriers (A = 1, r\ = 1/2) for k B T = W~ 2 E«, E p = 25E tT , 
g — 0.9, g p — 0.8, g a = 1.15 corresponding to e CT /e = 0.99, 
e p/ e o — 1-42 and 2E a /eQ = 0.86. Right: gray scale. 



The condition for opening the transition \n, s) <-> 
| n + l,s ± 1), without involving CDW and SDW, is 

-r]eV < -E p h g ± E a s < (1 - r))eV . (41) 

with n g = n g — n g cs (n, 0) and 

Cr + C s /2 
Cl + Cr + C g 

For symmetric barriers rj = 1/2 and A = G^/G^ = 1. 
The scheme of the transitions consistent with (|41f) is 
shown in Fig. States with spin s and —s are degen- 
erate (cf. (jSJ). First, we consider transitions between 
the ground states, |n, 0) <-> \n + 1, 1) associated with the 
lines Eph g = (rj — l)eV and E p h g — rjeV, respectively. 
These lines divide the plane (V > 0, n ) in three re- 
gions. The shaded regions labeled n and n+ 1 denote the 
Coulomb blockade regimes with zero conductance. In the 
other regions, both of the ground states have finite occu- 
pation probabilities and the quantum dot is conducting. 

Since the state \n+ 1, 1) is now occupied, the transi- 
tion \n+ 1,1) — > | n,2) becomes available at sufficiently 
high voltages. The corresponding activation threshold is 
given by E p h g = rjeV - 2E a (cf. gTJ with s = 2). The 
transition channel |n, 2) — > |n + l, 1) is always open in 
the transport domain. Inside the domain where \n, 2) 
is occupied we can achieve the transition for |n, 2) — > 
|n+ 1,3) by increasing V. By iterating the procedure 
with increasingly higher values of the spin we get the 
fish-bone pattern shown in Fig. 

So far we only have considered the excited states with 
higher spin of the dot. In order to complete the picture, 
we must include all the transitions involving collective 
charge and spin excitations. This enhances considerably 
the complexity of the spectrum, since at high enough 
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FIG. 9: Scheme of transition lines in the differential con- 
ductance for symmetric barriers (rj = 1/2, A = 1). Shaded 
regions denote Coulomb blockade. Thick black lines: ground- 
state to ground-state transitions; dark-gray sequence of lines: 
transitions to higher-spin states (spin values of the states in- 
volved in the transition are indicated); dashed lines: SDW 
and CDW excited states. For simplicity we have dropped 
indices n, n+1 denoting the transitions by spin values s only. 




FIG. 10: The differential conductance G (arbitrary units) as 
a function of V (units E^/e) and n g for asymmetric barriers 
with A = 50, C L = 5C R , C g /C R = 0.01, k B T = 10" 2 £ CT , 
E p = 25E a , go = 0.9, g p = 0.8, g a = 1.1 corresponding to 
e CT /e = 0.995, e p /e = 1.37, 2E a /e = 0.9. Several transi- 
tions associated with NDC (white lines) parallel to transition 
lines |n+l, 2s— 1) — > \n, 2s, p, a) are observed. Right: gray 
scale. 



bias voltage each transition of the type \n, s) — > \n',s'} 
can also occur via the channels \n,s) — > \n' , s' ,l p ,l a ) 
\n',s'). For example, consider the CDW- and SDW- 
channels for the transition |n, 0) — > \n + 1, 1). These cor- 
respond to the equations E p h g = (77 — 1) eV + l p e p + l a e cr . 
The lines are parallel to the line \n, 0) — * \n + 1, 1). Anal- 
ogously, including CDW- and SDW-channels in the tran- 
sition \n + 1, 1) — > |n,0) gives rise to E p n g — rjeV + 
l p £ p + lo£ G and to transition lines parallel to the line 
I n + 1,1) -» |n,0) (Fig. 13. 

B. Non symmetric barriers 

Figure lTUI shows the differential conductance for asym- 
metry A = 50. 

The main difference between the results in Fig. 1101 and 
those obtained for symmetric barriers (cf. Fig. [SJ) are the 
negative conductance features associated with the white 
transition lines. They are due to the presence of inter- 
action, e p 7^ E a , as can be seen from Fig. 1111 where the 
differential conductance is shown for g p = g a = 1. 

For V > electrons flow from right to left. Then, for 
A > 1 the electrons traverse a higher barrier tunneling 
into the dot, and a lower barrier tunneling out of the dot. 
Thus, states \n, s) will have a higher occupation probabil- 
ity as compared to states \n+ 1, s + 1). This "trapping" 
phenomenon at sufficiently large asymmetries can cre- 
ate a bottleneck for the electron transport and enhances 
the probability to have a NDC. The occupation proba- 
bilities for the lower spin states are shown in Fig. 1121 for 
increasing bias voltage. If the occupation of a higher-spin 
state e.g. (n, 2) is favored on the expense of the ground 
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FIG. 11: The differential conductance G (arbitrary units) as 
a function of V (units E a /e) and n g for A = 50, Cl/Cr = 
5, C g /C R = 0.01, k B T = 10- 2 E a , E p = 25E a , g = 0.9, 
dp — 9<y — 1-0, 2E(j = e CT — e p = so- Despite the strong 
asymmetry there is no NDC. Lines parallel to the transitions 
\n + 1, 1) — > |?2,0} tend to have very small intensity. Gray 
scale as in Fig. 1101 



state, (n, 0) , it accumulates occupation probabality and 
this eventually leads to a reduction of the current. 

However, it is important to note that trapping alone 
is not sufficient to induce NDC. If the system does not 
exhibit spin-charge separation (Fig. Ill|) no NDC occurs, 
independent of the strength of the asymmetry. To un- 
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v 4 

FIG. 12: Occupation probabilities P„ t o (dashed), P n ,2 + Pn,-2 
(dashed-dotted), P n +i,i + Pn+1,-1 (dotted), and current I 
(full line, units eT^ ) as a function of the bias voltage V 
(units E a je) for n g = 0.457; other parameters as in Fig. 1101 



derstand this point better, we use a simple model that 
can be treated analytically. We assume g p = g„ = go = 1 
and fceT <C e a . 




FIG. 13: Transport without spin-charge separation. Left: 
transition lines for asymmetric barriers (A > 1) and g p — 

g a — 1. Right: values for the current / in units To = sTq 
for A — > oo and T — 0, keeping Gr finite. 



We refer to Fig. [23 (left) where several regions of given 
higher-spin states are shown: {^4 m } with s max = 1, {B m } 
with s max = 2. In these regions the transition rates occur 
in a well defined pattern. For instance, for {A m } we have 

r^ 1 = mrf and = T h) , while for {B m } we have 

r& = t?\ r< L ] = (m + i)if \ r& = mr^ and 

= 2Fq . This allows to solve the master equa- 
tion and evaluate the currents (cf. Fig. ^3 right) for 
A — > oo. This limit is performed by assuming Gr, fi- 
nite and corresponds to the most favorable situation for 
generating NDC. One finds that the conductance lines 
parallel to |n+l,l) — > \n, 0) vanish. The lines paral- 
lel to |n, 0) — > |n + l, 1) correspond to PDC. Thus, for 
A —>■ oo, the absence of charge-spin separation leads to 
a landscape of PDC peaks. The numerical results for 
A = 50 (Fig. Ill|) are consistent with this picture. We 
have confirmed numerically that the above results remain 
valid also for interacting leads, go < 1, and for T ^ 0. 



C. The five-states model 

Among the transitions in Fig. ^3 which show NDC, 
we will concentrate on those in the gray region for the 
transition \n, 2) — > |n + l, 1) (cf. Fig.fTl|l. This is divided 
into three regions I, II, III depending on the presence of 
SDW and CDW states, l a = l p = 0; l a = 1, l p = 0; 
l a = 0, l p = 1, respectively. Here, only the five states 
|n,0), |n + l,±l), |n,±2) contribute to transport. For 




FIG. 14: Regions analyzed analytically with respect to the 
NDC: I, II, III correspond to l a = l p = 0; l a = 1, l P = 0; 
l a = 0, l p = 1, respectively, in the transition |n, 2) — > |n+l, 1). 

these, the master equation is solved analytically. The 
current is 





p(L) r (L) " 
1 1^0 1 1^2 


F (R) 
1 2-»l 


r (L) r (R) , 9p (R) 
1 l-i-O 1 2^1 ' zl 0^1 


r (R) p (L) " 
1 2^1 1 1^2 



(42) 



with the transition energies AUq^i = eV — 2E a , 
AU 2 ^i = eV, AUi^o = 2E a and AUi^ 2 = inside the 
corresponding rates. For k B T <C E a and go not too small, 
one can assume that along the lines \n + 1, 1) — > |n,2), 
|n +1,1) -> |n, 0, 0, 1), |n + 1, 1) ->■ |n, 0, 1, 0), and not 
too close to the crossing points, the only rates that con- 
tribute to the derivative with respect to V are and 
r^o- With this the differential conductance is 

G = W £ KdvT [ H 1+p {V) , (43) 
p=±1 

where </> = 2r[ ) I 2 i r2-i 1 and 

a p = ^+p(rW_ 2r w i )rB 1 _ p 

2>o = r( L 2 r2 R J 1 + 2r^ 1 (rS 1 + ri L 2 2 ) . (44) 

Expression i|43|) can change sign depending on the fac- 
tor A p . For G along |n + l, 1) — > |n, 2,l p ,l a ) we have 
p = 1 , and p = — 1 for G along |n + 1, 1) — > |n, 0, Z p , /„). 
Consider the transition \n + 1, 1) — » |n, 2), thus p = +1. 
Along this line, has contributions from the ex- 

cited states, while the other rates are Fq L ' R ' ) . For non- 
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interacting leads, go 
lent to 



1, the condition Ai < is equiva- 



r, 



(R) 



2—1-1 



A> 1 



We define A c as the critical asymmetry above which 
NDC is found. In region I, = 2. In region II, A^ = 
2 (2g a + 1) / {2g a - 1) > A§> . Thus, for A > A< n) one 
has NDC in both regions I and II. On the other hand, 
the condition for having NDC in region III is 



A ( 2g g a - g 



9 a 



2 \ 2 9 P 9a +9 P + 9o 



> 1. 



(45) 



For small exchange interactions this requirement cannot 
be fulfilled. Indeed for V cx < V /2 one has g 2 . + 3g 2 — 
1 " 2 g 2 > which cannot be satisfied simultaneously with 



We conclude that for sufficiently large asymmetry 
NDC is generally found in regions I and II and PDC in 
region III at the transition \n + 1, 1) — > |n, 2). 

The conditions for NDC along the lines |n + 1, 1) — > 
\n, 0, 0, 1) and \n + 1, 1) -> \n, 0, 1, 0) is 



1 



r (R) 
L _o 



A> 1 



It is clear that both in I and in II no NDC are found. 

(R) (R1 

In zone III, however, it is possible to have > 2T 

with NDC for 

A > = 2 -9p + ja + 2g P g CT 



9 P 



^9 P 9a 



If A > max^^A^,^ 11 ^} the NDC-PDC pattern 
has therefore a "photographic- negative" -like shape. This 
means that if the line |n + 1, 1) — > \n,2) (p = +1) has 
NDC, the two adjacent ones (p = —1) have PDC, and 
vice versa. This feature can be found in many regions of 
the density plot for not too high bias voltage (cf. Fig. [7J • 
These results indicate that in addition to the condition 
£p 7^= £<r one must have non-integer steps in the transition 
rates as a function of the energy in order to find NDC. 
It seems that this is a genuine feature in present case of 
an interacting system with non- Fermi liquid correlations, 
i.e. with spin-charge separation. 



D. Spin-flip relaxation 



with s' = s ± 2, n' — n, and consider w to be an esti- 
mate of the typical spin-flip rate. This spin-flip mecha- 
nism acts against the trapping of the higher spin states, 
providing an escape possibility \n, 2) — ► |n, 0). This is 
expected to decrease the chance to have NDC. 

The differential conductance for different values of the 
spin relaxation rate w is shown in Fig. 1151 for voltages 



within the five-states region, s n 



±2. The peaks for 



a positive bias are reduced by increasing w, especially 
the NDC state which eventually becomes positive. For 
negative bias, electrons traverse a lower barrier tunnel- 
ing into the dot, and a higher barrier tunneling outside: 
Pn,a,Pn,±2 < Pn+i.±i- Therefore the relaxation from the 
states |n, ±2) is almost ininfluent, since these states have 
already a small occupation probability. 




FIG. 15: Differential conductance (units e 2 Y^ /E a ) a function 
of positive and negative voltages V (units E a /e) in the five- 
states region for n s — 0.485 for spin relaxation rates w = 0.2 
(solid), 0.3 (dashed), 0.5 (dotted) and 0.6 (dash-dotted) in 
for k B T = lO" 2 ^, E p = 25£ CT , A = 10, 
0.9, g p = 0.7 and g a = 1.25. 



1/2, go 



units T^'(E„/u!. 



n 



As an example, we consider the above five-states model 
including spin-flip relaxation. The conductance has the 
same form as in (|43f) 



G =^ E ^d v T^_l 1+p (V) 



(47) 



but now 



p=±i 
i(R) r„„ , r( R ) 



2rni(™ + rni) and 



Since the NDC is related to the occupation of states 
with spins higher than that of the ground state, it is 
of great interest to introduce into the master equation 
an extra spin-flip relaxation rate r™_^ s ,. This will give 
insight into the depletion behavior of the NDC when spin- 
flip processes are present. We assume 



w exp 



< \s\ 
> \s\ 



(46) 



A w = 

-lip 



p(R) _ np(R) 



r (L) 
1 



r (L) r (R) 
1 1^0 A 2->l 



r (L) 



2r 



(R) 
0->l 



,(R) 



,(L) 



r (L) 
1 l->2 



2r 



(R) 



(48) 



As above, the sign of the conductance is solely given by 
Ap which is now a function both of the asymmetry and of 
the relaxation w. We concentrate on the region I (2E a < 
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eV < e a ) of the line |n + l, 1) — > |n, 2). The condition 
Af < reduces at T = to 



,(L) 



1 - 



1 



eV 



1 / eV 



2 A \2E, 

with a given in 127|) . The critical value u> c for the 
0.6 




FIG. 16: Phase diagram for the critical value of the spin-flip 
relaxation rate w (units (E a /Lo c ) a ) in region I as a func- 
tion of the asymmetry A for T — and different interaction 
parameters in the leads: go — 1.0 (solid), go = 0.9 (dashed), 
go — 0.75 (dotted) and g = 0.65 (dash-dotted). 



crossover is plotted in Fig. EH as a function of A. At 
a given asymmetry, for w > w c the conductance peak 
crosses over to positive values. The stronger the inter- 
action in the leads the smaller is the threshold in the 
relaxation for destroying the NDC peak. 



VII. CONCLUSIONS 

We have investigated the non-linear current voltage 
characteristic of a ID quantum dot described by a Lut- 
tinger liquid of a finite length connected via tunnel bar- 
riers to interacting quantum wires. The system contains 
four characteristic energy scales: charge and spin addi- 
tion energies, and the energies needed to excite charge 
and spin density collective modes in the quantum dot. 
We have discussed in particular the effects related to the 
presence of the electron spin. 

In the linear regime, we have reproduced the parity ef- 
fect in the distances between the Coulomb blockade con- 
ductance peaks. In contrast to the case without interac- 
tion discussed earlier^ the temperature behavior of the 
peak positions shows signatures of the non- Fermi liquid 
nature of the leads. At temperatures higher that the spin 
addition energy, the peaks become equidistant. 

In non- linear transport, we have found NDC features 
that are related to the presence of spin-charge separation. 
They are connected to states with higher spins that par- 
ticipate in the electron transport. 

This effect is different from the spin blockade phe- 
nomenon discussed previously^ In the spin blockade of 
type I, a "trapping" argument is used to explain the ef- 
fect. However, in this case the only state that leads to 
trapping (and thus can decrease the current) is a state 



with the highest total spin and highest energy for a given 
electron number. Thus, only one NDC peak can be 
present in a given n n + 1 transition. The "trapped" 
status is achieved because of the Clebsch-Gordan coeffi- 
cients that are introduced in the calculation of the tran- 
sition rates. 

In contrast, in the present case, each state with even 
(odd) spin can become trapped if (A — 1)V > (< 0) 
such that the transition lines that populate such a state 
exhibit NDC. Whether or not such a state is eventually 
trapped is determined by the transition rates which are 
here microscopically evaluated using the Luttinger liquid 
model. The result is that even if the system contains 
asymmetric barriers that lead to trapping of the appro- 
priate candidate states in the dot, it is still necessary to 
have charge-spin separation in order to obtain NDC. It is 
clear from the above that the effect of spin-charge sepa- 
ration is not only the trivial removal of energetic degen- 
eracy of the CDW and SDW states. The spectral weights 
contained in the tunneling rates l|25|l play a crucial role 
in determining whether or not a particular transition line 
corresponds to NDC. 

Since NDC is driven by asymmetry, there must be a 
threshold value A c above which the phenomenon occurs. 
This critical value is not universal. Each transition line 
has a different critical value for the asymmetry. 

We have assumed that the collective excitations that 
do not alter the z-component of the spin in the quan- 
tum dot have infinitely short relaxation times as com- 
pared with the excitations which are associated with spin 
changes. We feel that this assumption is justified in view 
of recent results suggesting that in quantum dots the re- 
laxation of states without flipping spins can be orders of 
magnitudes shorter than that associated with spin flips. 41 
Nevertheless, the question whether or not without this 
assumption the predicted NDC-phenomena would dis- 
appear is legitimate. In order to answer this question 
quantitatively, it is necessary to repeat the calculations 
including the collective states of the quantum dot as dy- 
namical variables in the master equation^ 

Spin-flip relaxation processes only seem to weaken ex- 
isting NDCs, as indicated by the results in the section 
IVIDI Thus, we can expect that by removing relaxation 
processes and including additional stable states in the 
transport process without flipping the total spin will not 
change qualitatively the trapping mechanism described 
in the section lYl Bl The NDC predicted here, which ap- 
pears to be a consequence of the trapping mechanism and 
the separation of energy scales for spin and charge exci- 
tations together with the Luttinger liquid features that 
enter the rates, will not be depleted. 

In experiment, several possibilities for measuring the 
predicted non-linear phenomena exist. The necessary 
condition for applying the above model is 



k B T < e < E F , 
with £o — -KVp/a- This is equivalent to 
At > a » Af 



(49) 



(50) 
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where At = ttwf/^b? 1 is the thermal length and Ap = 
27r/fcp the Fermi wave length. Tunnel barriers will in 
general be asymmetric in any case. 

In semiconductor based quantum wires, quantum dots 
are prepared via depicting the electron density such that 
eventually an electronic island is formed accidentally be- 
tween two impurities. 20 In this system, the Fermi energy 
appears to be relatively low, Ep sa 2meV, due to the 
almost depletion of the lowest sub-bandiSi The dot en- 
ergy level spacing is about 1 meV for a quantum dot of 
length a 0.2 /jm and an effective mass m* rs 0.07 m G 
(for GaAs). These parameters at the first glance seem to 
be outside the above region of validity of the model, al- 
though the condition with respect to the temperature 
is easily fulfilled. For a quantum dot of larger size 
(awl fim) the situation would be much more favor- 
able. This seems to be achievable in other semiconduc- 
tor based quantum dot systems, such as those fabricated 
by scratching techniques, 45 With this one could fabri- 
cate an appropriately scaled quantum dot in a quantum 
wire. Parameters to be achieved should be: width of wire 
« 50 nm in order to make a Fermi energy of a few meV 
achievable; distance between the tunnel barriers (point 
contacts) that define the quantum dot « 1 /urn giving 
Eq » 0.1 meV. Temperatures should be well below IK. 



With this technique, one would have the advantage of 
being able to adjust the asymmetry of the tunnel barri- 
ers. 

In carbon nanotubes the situation with respect to the 
energy scales seems to be even more favorable since Ep ~ 
2eV, ~ 8 • 10 5 m/s and the level spacing eo ~ 5meV 
with temperatures of the order T ss lOOmK. 4 However, 
here the non-interacting energy spectrum consists of four 
branches including the spin. The present theory has to be 
adjusted to this case in order to apply the results. Again, 
due to the fact that the predicted phenomenon seems to 
be quite generally valid, one can expect that experiments 
on carbon nanotubes should show NDC associated with 
the higher spin states. 
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